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Restrictions to molecular motion by barriers (membranes) are ubiquitous in biological tissues, porous 
media and composite materials. A major challenge is to characterize the microstructure of a material 
or an organism nondestructively using a bulk transport measurement. Here we demonstrate how the 
long-range structural correlations introduced by permeable membranes give rise to distinct features 
of transport. We consider Brownian motion restricted by randomly placed and oriented permeable 
membranes and focus on the disorder-averaged diffusion propagator using a scattering approach. 
The renormalization group solution reveals a scaling behavior of the diffusion coefficient for large 
times, with a characteristically slow inverse square root time dependence. The predicted time 
dependence of the diffusion coefficient agrees well with Monte Carlo simulations in two dimensions. 
Our results can be used to identify permeable membranes as restrictions to transport in disordered 
materials and in biological tissues, and to quantify their permeability and surface area. 



Brownian motion in a uniform medium is characterized 
by a single parameter, the diffusion coefficent D, which 
is a measure of mean square molecular displacement. A 
packet of random walkers spreads with time t according 
to a Gaussian distribution with variance (a? 2 ) = 2Dt. 
Complexity in the microscopic structure of a sample, such 
as heterogeneity in diffusive properties and restrictions to 
molecular motion, results in non-Gaussian evolution 1 ' 2 . 
In particular, the diffusion coefficient D(t) itself becomes 
time dependent. 

A fundamental question is how specific complexity fea- 
tures manifest themselves in the coarse-grained disper- 
sive dynamics. Random drifts are known to drastically 
slow the dynamics down to (a; 2 ) ~ ln 4 i, thereby sup- 
pressing the diffusion, D\t-+oo — 0, in one dimension 3 ' 4 . 
Their effect in higher dimensions is less profound 5-9 . 
Randomness in local diffusion coefficient preserves Gaus- 
sian diffusion at t — > oo, resulting in the finite limit 
Doc = D(t)\ t=00 . However, it causes a power-law disper- 
sion ~ ui d / 2 in the real part of the velocity autocorrelation 
function V{uj) — J£°dt e lLJt (v(i)i>(0)); this power law de- 
pends on the spatial dimensionality d and results in long 
time tails in the system's current-density response kernel 
V(t). w Studying dispersive diffusion, therefore, is a way 
to characterize the type of disorder and of restrictions to 
molecular motion in a complex sample. 

Here we consider the important class of restrictions to 
diffusion, the permeable barriers (membranes). Practi- 
cally, membranes play an essential role for the transport 
of ions 11 , water molecules 12-1 5 , and gases 16 in biological 
tissues 17 , and in porous and composite materials 18-23 . 
On a fundamental level, while occupying vanishing vol- 
ume, permeable membranes introduce important long- 
range correlations into the structure of a disordered sam- 
ple. We show that these correlations give rise to distinct 
transport features, qualitatively different from those in 
"uncorrelated" disordered systems 10 . 

We develop a minimal model which appeals to any d- 
dimensional complex medium where the restrictions by 



permeable membranes play a dominant role. We assume 
that diffusion is restricted by randomly placed and ori- 
ented infinite flat membranes, such as in the d = 2 exam- 
ple of Fig. 1. We find the diffusion coefficient as a func- 
tion of time or frequency, for all membrane concentra- 
tions and permeabilities. We demonstrate that the spa- 
tially correlated disorder introduced by the membranes 
results in a long-term memory that manifests itself in a 
distinct non-analytic low-frequency dispersion of the dif- 
fusion, T>(oS) — Dqo ~ y/u. The latter is equivalent to 
a characteristically slow decrease of the diffusion coeffi- 
cient, D(t) — Doo ~ 1/v? as t — >■ oo, Fig. 2, causing 
the mean square molecular displacement to increase as 
(a: 2 ) ~ 2Z?oot + const • yft. We relate this power law 
dispersion to the anomalously strong fluctuations of the 
amount of restrictions in a given volume, caused by the 
spatially extended nature of the disorder. This makes the 
yfbj ~ 1/ yi dependence present in any spatial dimension 
d as long as the membranes are flat on the scale of the 
diffusion length, as confirmed numerically for d = 2. 

Our finding emphasizes the role of spatial correlations 
in disordered samples in contrast to the short-range dis- 
order typically considered in classical 10 and quantum 24 
transport. It constitutes a novel disorder universality 
class of classical random media. The y/uj ~ 1/Vt depen- 
dence in the system's low-frequency dynamics can serve 
as a unique "fingerprint" of the permeable membranes 
within the complexity of realistic samples in which the 
transport can be studied via the electrical 25 or heat con- 
duction, or by the diffusion- weighted NMR 14 . 



Model 

We begin by outlining a minimal model of a sample in 
which the dominant restrictions to molecular motion are 
random permeable membranes. A membrane is an ide- 
alization of a thin slice of a poorly diffusive material, as 
long as its thickness is negligible compared to the short- 
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est observable diffusion length. In the limit when both 
the diffusion coefficient D m of the membrane material 
and its thickness l m vanish, the ratio k = D m /l m is by 
definition the permeability. The effect of a membrane is 
described by the boundary condition 17 ' 26 ' 27 
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Here Do is the unrestricted diffusion coefficient. The con- 
dition (1) means that, at each point r m of the membrane, 
the density ip of random walkers experiences a jump pro- 
portional to the component of the current J along the 
normal n to the membrane surface. 

The permeability has the dimensions of velocity. In 
what follows, we find it useful to associate an effective 
thickness, 21 = Dq/k, with a membrane. This length 
scale is defined relative to the free diffusion coefficient 
Dq. Its physical meaning is derived from the condition 
(1) and is illustrated in Fig. 3a; the membrane indeed 
appears DQ/D m ^> 1 times thicker than its "nominal" 
vanishing thickness l m . 

Consider now diffusion in a macroscopic sample em- 
bedded with multiple randomly placed membranes, each 
one imposing the condition (1). The diffusion propaga- 
tor depends on their number, shape, and spatial distri- 
bution. The number of membranes is characterized by 
the ratio S/V of their total surface area to the sample 
volume. Here we adopt the convention from the porous 
media literature 15 ' 20 : a membrane has two faces, so that 
the membrane's surface area is counted twice in S. The 
shape and the spatial distribution of the membranes vary 
greatly depending on the physical context. However, as 
it will follow from our treatment, the main dispersive fea- 
tures of transport can be captured by making the sim- 
plest assumption which also allows us to keep the number 
of parameters to a minimum. Namely, below we consider 
the membranes as infinite d— 1 dimensional planes placed 
and oriented in a completely random (uncorrelated) way, 
dividing the sample into pores with random shapes, as 
shown in Fig. 1 for d — 2. In this case, the ratio S/V is 
all what is needed to characterize the geometry. 

In what follows, we first outline our main results for 
the model medium. Next, we derive them, compare with 
numerical simulations, discuss and generalize. 



Results for randomly oriented flat membranes 

In the t — ¥ oo limit, the diffusion becomes Gaussian 
with the reduced diffusion coefficient 
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Here the dimensionless parameter £ quantifies the ability 
of membranes to hinder the diffusion. It is the "volume 
fraction" occupied by the membranes based on the effec- 
tive thickness I defined above. The result for is exact 
in d = 1, ref. 28. As follows from Fig. 2, it is a good ap- 
proximation for d > 1, with 1/d = ^cos 2 9} in equation 




FIG. 1: A fragment of a two-dimensional patch with ran- 
domly placed and oriented membranes (one of the disorder 
realizations used in the simulations). 



(2) arising from the mean-field orientational averaging, 
cf. Methods section. Strong restrictions correspond to 
£ > 1, when the domains of thickness ~ £ associated 
with each membrane overlap, such that transport across 
a membrane is affected by its neighbors. 

For finite t, we shift to the frequency representation 
since the restrictions are stationary. Technically, we focus 
on the diffusion propagator G averaged over the disorder 
in positions and orientations of the random membranes. 
After disorder averaging, the propagator becomes trans- 
lation invariant. Its pole in the frequency - wave vector 
representation, G~q = — iuj +T>(uj)q 2 + O (q 4 ) , defines the 
dispersive diffusion coefficient T>(u)) which is a retarded 
response function that relates the disorder-averaged par- 
ticle current J u r = —T>(u>)d r ip UJir to the density gradient. 
The corresponding time-dependent diffusion coefficient 
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is given by a contour integration along the real axis with 
all the singularities in the lower half-plane of complex uj. 

We find T>(uj) in three steps. First, we develop a scat- 
tering approach for the transmission events described by 
the boundary conditions (1) at each membrane. Sec- 
ond, we solve the problem perturbatively in the volume 
fraction £ at the mean-field level, starting from infinitely 
permeable limit T>(lo) = D Q . These two steps are done in 
the Methods section. The dispersive diffusion coefficient, 
valid up to 0(C), is 



Here z<. 



T>(u) = D 1 



and 



c 



1—ZtL 



c«i. 



(4) 



e/D = D /(2 K f 



3 





2.5 



8 



1.5 



a 1 



0.5 



O O O O o O 



10 



10 



C 



10 



FIG. 2: Time-dependent diffusion coefficient D(i) for the two- 
dimensional random medium of Fig. 1. a, Comparison of 
the RG solution (5), red, with the Monte Carlo simulations, 
blue, for the set of decreasing permeabilities, corresponding 
to C = 0.5, 1, 2, 4, 10, 20, 40, 100, 200, 400 (top to bottom). 
The diffusion time td is marked by blue circles, b, Scaling be- 
havior of D(t). As the strength £ of the restrictions increases 
from bottom to top, the numerical curves begin to collapse 
as a signature of the universal behavior (10). Dashed lines 
show the £ — > oo limits from our RG solution: the "imperme- 
able" limit D(t)/Do = 2-TD/i and the scaling limit (10) with 
62(00) = 4/- v /7r. c, Parameters of the scaling limit (10), C2(C) 
(filled circles) and D 00(C) (open circles), determined from the 
fit of the simulations in b to equation (10), compared with 
the RG predictions (solid and dashed red lines). Thin solid 
line is the RG limit 6*2(00) = i/y/n. 



At the third step we extend equation (4) onto the non- 
perturbative domain C ^ 1 using the real-space renor- 
malization group to account for the multiple transmission 
events self-consistently, as described below. Thereby we 
obtain our main result for all C, 
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Effective circuit 

Before deriving our main result (5), let us note that the 
perturbative limit (4) can be represented as a simple ef- 
fective circuit. As the random walkers are uncharged, the 
current J has only the diffusion component, and the dis- 
persive diffusivity T){uS) defines its response to a density 
gradient d x ip rather than to a potential bias. With this 
important distinction, consider interpreting T>(u>) as the 
analog of the conductivity, in the spirit of the Einstein 
relation, albeit at finite frequencies 25 . This way, the 0(C) 
result (4) can be viewed as an "impedance" Z(u>) of a hy- 
percube V = L d , L d - 2 Z(uj) = 1/V(uj) l/D a +nZ m (w), 
which acts as a one-dimensional disordered transmission 
line of length L with point impedances 



Z m (u) = , £fm(w) = k(1 - z u ) 
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placed at random positions in series, with density n — 
S/(2Vd). This addition of independent impedances, a 
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FIG. 3: a, Meaning of effective membrane thickness: 21 — 
Do/k is the distance by which one should shift the density 
profile ip(x) on each side of membrane if one were to heal 
the jump discontinuity (1). b, The equivalent circuit for 
membranes embedded in a medium with diffusivity Do. The 
membrane "impedance" can be represented as a resistor with 
conductance n in parallel with the permeability-independent 
Warburg element W with conductance gw{w) = \/iu)Do/2i. 
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consequence of the mean-field description (see Methods 
section), will be justified later, when we discuss the origin 
of the y/to dispersion. 

The equivalent circuit for a single membrane (6) is 
shown in Fig. 3b. In the dc limit, a membrane acts 
as a resistor with the "conductance" k = g m (0); as 
( = 2ln 7 each membrane effectively adds the length 21 
to the original clean wire L if its dc resistance were to 
match Z(0) oc At finite frequencies the resistor 

is shunted by the Warburg element 29 with the "conduc- 
tance" gw{w) = —kz u = ^iujDo/2i independent of the 
permeability. The Warburg impedance, first observed at 
a flat metal electrode in an electrolyte by Kohlrausch 30 
and Wien 31 in the 19th century, is associated with the 
diffusion- limited response 29 . Incidentally, the impedance 
Z(lo) is of the Cole form 32 ' 33 with power law exponent \. 

Renormalization group 



Integrating equation (8) yields our main result (5). 
In the dc limit X>(0) = D^, the RG flow 1/Doo|l' = 
1/DqoIl + (N' — N)/(nLd) recovers the above exact re- 
sult for d = 1, if one identifies L with the microscopic 
scale such that £>oo|l = D (the bare diffusivity) and 
(N' — N) /Ld — > n, the final membrane density. Written 
as RG flow, however, the evolution of the diffusivity is not 
tied to a particular microscopic cutoff; indeed, the diffu- 
sivity -Dqo|l can itself originate from "more microscopic" 
membranes on the scale finer than L, and so forth. 

We also note an interesting equivalent way of looking 
at the RG flow: Formally identifying the particle number 
density ip entering the diffusion equation with an electro- 
static potential, the dc limit corresponds to the electro- 
static screening of the one-dimensional electric field (the 
slope d x ip) by the dipoles (membranes acting as dou- 
ble layers). Coarse-graining corresponds to more dipoles 
contributing to the screening 35 , which renormalizes the 
effective dielectric constant, 1 — > 1 + £. 



To calculate the response T>(us) of the disordered trans- 
mission line for finite £, we employ the following scaling 
argument, which we first develop for d = 1. Consider a 
slab of length L with the diffusivity V{uj)\l. Let us ex- 
tend the slab to the length L' = bL, b > 1, and rescale; 
now the slab's length L' in new units is back to L, while 
the original, shorter slab has length L — > L/b. The 
rescaled slab conductance decreases due to adding extra 
N' — N = NSb membranes, 5b = b — 1. The diffusivity is 
then reduced as 
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according to equation (4), as long as the added mem- 
brane density Sn = (N/L)Sb is small, 5b <C 1. Since we 
assumed from the beginning that the membrane positions 
are uncorrelated, adding a small number of membranes 
in an uncorrelated way at each step is consistent. Choos- 
ing an infinitesimal dn oc db, we represent equation (7) 
in the differential form, obtaining the real-space renor- 
malization group (RG) equation 

dn " K-\yJiufD(u) ' 

This is a telegraph equation for a disordered transmis- 
sion line. Alternatively, it can be obtained by adding 
membranes in small increments Sn in a macroscopic 
sample of a fixed length with bare diffusivity Dq, and 
applying the relation (4) treating all previously added 
membranes in the effective-medium fashion at each RG 
step. Either way, the diffusivity flows to lower values 
as long as the added membranes have finite permeabil- 
ity k < oo. In d dimensions, rescaling a hypercube L d 
as in ref. 34 is equivalent to rescaling a d = 1 slab in 
which the effective one-dimensional membrane concen- 
tration n — S/ (2Vd) . At the mean-field level, the prob- 
lem is always one-dimensional, with the fraction 1/d of 
membranes exerting full resistance. 



Regimes for D(t) and numerical simulations 

Here we discuss the time dependence (3) of our result 
(5) and compare it in Fig. 2 with Monte Carlo simulations 
of diffusion in the two-dimensional geometry of Fig. 1. 

First we note that, for £ ~ 1, the agreement between 
theory and simulations is very good, which illustrates 
that the RG framework is able to extend the perturbative 
result (4) to the case of moderate disorder. Next, we turn 
to the more interesting and general case of the strong 
disorder, £ 3> 1, where there are three distinct regimes in 
D(t), separated by the two time scales, r r » T£>, defined 
below. 

The initial decrease of the diffusion coefficient 
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reproduces the well-known universal short-time \fl- 
expansion 20 . The correction to free diffusion arises from 
the random walkers within the diffusion length <~ y/Dot 
from the barriers; the term oc t is the permeability 
effect 23 for the case of flat pore walls. The limit (9) is 
valid as long as only a small fraction of random walkers 
has encountered the pore walls, i.e. for t <C t d , where 
td = a 2 /2Dq is the diffusion time across the typical pore 
size a ~ 1/n = 2dV/S. Equation (9) follows from the 
first two terms of the \z u \ ^ £ expansion of (5) and sub- 
sequent application of (3). This expansion can be per- 
formed using the perturbative limit (4), as it becomes 
universally valid for large frequencies, lutjj 3> 1, when 
only the walkers within the wavelength <~ ^/D /w from 
the membranes are restricted. There is a perfect agree- 
ment between theory and simulations in this regime for 
all C 

Highly restrictive membranes look completely imper- 
meable for t <S! T r , where r r = V/(nS) is the res- 
idence time 1 within a typical pore; in our notation, 
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r r = ((/d)T£>. As a result, the mean square displace- 
ment (x 2 ) is bounded by ~ a 2 , and D(t)/D <~ r D /t. 
For our RG result (5), the "impermeable" behavior oc- 
curs for 1 <C y/( -C \z u \ <C Cj and yields the purely imag- 
inary (-independent limit V(oj) ~ —iua 2 (similar to the 
purely reactive electrical conductivity) , corresponding to 
D(t)/D = 2r D /t, Fig. 2b. Both our result (5) and the 
simulations display the 1/t dependence for ( > 100, with 
the theory somewhat overestimating the D(t) obtained 
from the numerics for very large (. The discrepancy be- 
tween the RG result and the simulations is maximal for 
t <~ T r as expected from the discussion of our effective 
medium approach (see the Methods section). 

For times longer than r r , the system becomes aware 
of the finite leakage across the membranes, with D(t) 
approaching the limit (2) which can be qualitatively esti- 
mated as hopping on a d-dimensional lattice with a step 
a over the time t t , Doo ~ d 2 /(T r d). 

Remarkably, the way D(t) approaches slows down 
from ~ 1/t to 1/v^-decrease in any dimensionality d, 
defining the novel disorder universality class, represented 
by the randomly placed infinite membranes: 



£>«, \l + C d (C)\lj 
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This follows from V{uj) = D^ + 0(y/uj) in the limit 
\zw\ "C The regime (10) can be clearly seen from 

the simulation results of Fig. 2b where time is in the 
units of r r for a set of permeabilities. In this limit, the 
relative deviation [D(t) — D^/D^ < 1. We empha- 
size that for £ — > oo, the dependence (10) approaches a 
universal scaling law with a fixed d = C<<(oo), which 
can be represented only in terms of the effective parame- 
ters .Doo and r r , independent of the original microscopic 
parameters D , n and S/V. Approaching this law corre- 
sponds to the collapse of the simulation curves in Fig. 2b 
onto one universal curve for £ > 100 and t > r r . Our 
RG approximation, yielding C'd(oo) = yj&d/ir, overesti- 
mates the 6*2(00) obtained from the fit of the numerical 
curves to equation (10), as illustrated in Fig. 2c. For 
finite C, the parameters given by equation (2) and 
C d (0 = C d (oo) VC (\/r+C - 1) /(l + C) agree with the 
simulations within about 15% and 30% correspondingly, 
even in the strongly non-perturbative regime of £ ^ 100. 

The fact that the RG solution works fairly well even 
when the "small parameter" ( > 1 can be understood 
in terms of the flow to the moderately disordered limit 
with strongly renormalized parameters, D n — > and 
t — > r r . Indeed, for t >• r r , the medium effectively 
looks as if it had a decreased diffusion coefficient -Dqo. 
This, in turn, reduces the contrast between the diffu- 
sion coefficient inside the membrane and of its surround- 
ings, effectively changing Do — > D^ in the boundary 
condition (1), thereby reducing its effective thickness 
i — > l = _D 00 /(2k) ~ £/( ~ a down to the mean distance 
between membranes. Hence, the renormalized disorder 
strength is reduced, ( — > ( = SijVd ~ 1. Likewise, the 



renormalized time scale t — > f = £ 2 /Dq ~ r r matches 
the residence time. At this point the perturbative limit 
(4) is matched. 

Conversely, when the membranes are very permeable, 
( < 1, the residence time in a pore is of the order of td, 
and the time scale t t defined above becomes obsolete. 
In the perturbative limit £ <g; 1 our result (4) becomes 
exact, the intermediate 1/t regime does not appear, and 
the short-time behavior (9) directly crosses over to the 
nonanalytic dispersion of the form (10), 
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The role of r r is now played by the time scale t = 
((d/2)r r defined after equation (4). Hence it is now the 
time scale r -C r r which is associated only with the mem- 
brane properties and not with their density S/V, that de- 
termines the dispersion of diffusion for randomly placed 
highly permeable membranes, as discussed below. 



Origin of the y/Uj ~ l/\/t dispersion 

The defining feature of our model, the non-analytic 
dispersion T>{ui) — D^ ~ y/cj, is present already at the 
level of the perturbative result (4). We now discuss its 
physical origins, as well as justify the mean-field approach 
utilized to obtain equation (4) and its nonperturbative 
counterpart (5). 

In the perturbative £ <C 1 limit, the two relevant time 
scales are td and r = \C, 2 td- We first focus on the 
limit t <C td, where adding independent contributions 
(6), cf. equation (20), is valid since multiple scatterings 
off different membranes are not yet important. In this 
case, the dispersive behavior (11) for t <C t <C td is 
a result of scattering off individual membranes, in the 
following sense. Consider a constant density bias *o 
introduced across a membrane at t = 0. The single- 
membrane conductance (6) yields the transient current 
response J(t) = K x I / o(l — -Psurv(O) which involves the 
survival probability P SU rv|t»r = \fr~pKi of a random 
walker released a distance t = \/Dqt from the origin, 
to never return to the origin during time t. The frac- 
tion P SU rv(£) represents the walkers who have managed 
to wander around the "biased" side and have not yet en- 
countered the membrane. Hence, the transmission pro- 
cess is represented in terms of multiple attempts to re- 
turn to the position of the membrane with net prob- 
ability 1 — -P S urv(i), times the rate oc n to overcome 
the restriction. Equivalently, the conductance kernel 
9m{t) = J^ e- iut g m (uj) = np(t) is the Levy-flight prob- 
ability distribution function p(t) = — P S urv — ci~ 3 / 2 for 
time intervals t between successive returns of a Brownian 
particle to the or igin 2 , with r being the short-time cutoff, 
and c = ^/t/47t. Note that the contribution KP surv (t) of 
the "survivors" , originating from the Warburg contribu- 
tion gw to (6), is independent of the permeability, as 
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they have not yet encountered the membrane. In the dc 
limit t = oo the membrane acts as a simple resistor since 
in one dimension all Brownian particles eventually return 
to the origin. Incomplete transmission for r <C t <C td 
produces the \j\ft correction to the dc "conductance". 

For longer times t ^> td, multi- membrane scatter- 
ings determine the time dependence of D(t) which be- 
comes sensitive to correlations between positions of mem- 
branes in space. Indeed, at large diffusion times, the 
medium is effect ively coarse-grained into domains of the 
size <~ L(t) = \/1Dt >• a. These domains have slightly 
different diffusion coefficients Dj due to a different num- 
ber Nj of membranes falling into them. The diffusion 
coefficient D(t) = (Dj) is the ensemble average over the 
domains; in one dimension, its relation to is deter- 
mined by adding the "dc resistances", l/D^ = (1/Dj), 
leading to D(t) ~ D x + ({5 D) 2 ) / D x . The variance 
((SD) 2 ) ~ ((SN) 2 )/L 2 ~ 1/L(t) scales as i" 1 / 2 as a re- 
sult of the Poissonian statistics in the membranes' po- 
sitions. This simple argument elucidates the meaning 
of the previously obtained one-dimensional results 36-39 . 
(The crucial role of the Poissonian statistics is obvious 
from a qualitatively different 1/t decrease of D(t) in the 
case of periodic membranes in d = l. 27 ) 

Remarkably, our results (4) and (5) are compatible 
with the statistics of restrictions provided by randomly 
placed membranes in all dimensions d. The scaling 
((SD) 2 ) ~ ]-/L(t) persists for d > 1 as the random in- 
finitely long membranes generate qualitatively stronger 
fluctuations in the distribution of {Dj}, exemplifying the 
role of the spatially correlated disorder introduced here. 
Furthermore, the prefactors of the <~ y/co terms of 

equations (11) and (4) appear to be exact in d = 1, as it 
can be checked by generalizing the lattice calculation of 
ref. 36 onto our case of the membranes having random 
positions on a line, validating our results for 4> td- This 
is why the numerical simulations for small Q agree with 
theory very well even for very long t (Fig. 2). This ob- 
servation allows us to conclude that adding independent 
membrane contributions in a mean-field way is compat- 
ible with the Poissonian statistics of the disorder for all 
w. Formally this allows us to extend the above connec- 
tion with the return-to-origin probability onto t >• tu 
by substituting the rest of the system via a "featureless" 
effective medium with diffusivity D^. The equivalence 
between Poissonian disorder in membranes' positions and 
the return-to-origin probability persists in the RG solu- 
tion, as it is utilized at every RG step (7). As a result, 
these two pictures remain equivalent in the £ ^> 1 limit 
after the effective- medium substitution r — > r r , and the 
scaling behavior (10) is justified for all ( and t. 



Outlook 

In this work we have introduced a novel class of disor- 
der, represented by straight permeable membranes which 
are randomly placed and oriented. A random medium of 



this type has dispersive diffusion T>(uj) — D x <~ cj 1 / 2 in 
any dimensionality d. This dispersion for d > 1 is notably 
more singular than its well studied counterparts giving 
uj d l 2 in spatially uncorrelated random media 10, 36-39 . As 
the w 1 / 2 singularity originates from long-range spatial 
correlations introduced by the membranes, the disper- 
sive behavior of the form (10) will persist in a vari- 
ety of samples with random locally flat restrictions, and 
will become increasingly important for longer diffusion 
times when the contributions from shorter-ranged disor- 
der types vanish. However, the behavior (10) will eventu- 
ally terminate when the membranes cease to be flat past 
their intrinsic correlation (gyration) radius r c , i.e. for 
t > t c ~ r 2 /D(t c ), when the structural memory is for- 
gotten, and the dimensionality-specific dispersion ~ bj d l 2 
of a coarse-grained random medium sets in. 



Methods 

A single membrane at the origin. Consider first 
a one-dimensional problem. The boundary condition (1) 
can be represented as a local scattering term in the dif- 
fusion equation 

M = D d 2 x iP + V x iJ> , V x %l> = -2D £6'(x) [^VUo ■ 

(12) 

Here S'(x) is the derivative of the Dirac delta- function at 
the position of the membrane. 

We focus on the Green's function Gp, x ,x' of the prob- 
lem (12), defined by the initial condition Qt- x ,x'\t=o = 
5 (x — x'). As the problem is stationary, we shift to the 
frequency domain, dt — > —ioj. In this representation, the 
Green's function Gu-x,x' is formally an operator inverse 



(13) 



where the coordinates x play a role of indices in infinite- 
dimensional matrices, and the free propagator 



1 



-iuj + D q 2 ' 
The result of inversion (13) is the Born series 



(14) 



a ,-r (0) i / A™ r (0) v r {0) 

— ^ui-x-x' "r" / ux l ^u-x-xi v w;ii u u;ii-i' \ LO ) 

where the key quantity, the so-called "full vertex" satis- 
fies the Dyson equation, which in the Fourier representa- 
tion reads 

V^.kM = V kM + J ^ V k , q GMv u . qik , . (16) 

In our case, the bare vertex is separable, 

V^V ktk , =2D £kk'. (17) 



7 



The solution of equation (16), the full vertex 

2D £kk' 

Vuj;k.k> = 7== ■ (lo) 

1 — l\JlLOT 

Throughout this work, u) is understood as having an in- 
finitesimal positive imaginary part +i0 due to causality; 
this way the retarded response functions are analytic in 
the upper- half-plane of the complex variable lo. 

The w 1 / 2 singularity parallels that of the energy- 
dependent quantum scattering off a localized potential 
in one dimension, identifying iuj — > e + iQ. The differ- 
ence, however, is that in the quantum problem, energy 
comes with a negative power e -1 / 2 in the denominator of 
the full vertex: Scattering is maximum for high-energy 
particles, whereas low-energy ones do not notice a lo- 
calized potential. In our case the situation is opposite: 
A membrane is most important in the dc limit, with 
V w; fc ; fc' — > Vk,k', whereas at short times its scope is rc- 
duce d dow n to the particles within the diffusion length 
~ y/D /uj, with V u; k,k' ~ Vk^k'/V^J- 

Multiple membranes. Averaging over positions of 
randomly placed membranes with concentration n is done 
in a standard way 1 ' 2 ' 24 , by introducing the self-energy 
part Y^ W q, 

GZ] q = G£» _1 - . (19) 

On the mean-field level, we keep in only the exact 
interaction vertex with a single membrane, 

~ nV w;M . (20) 

The mean-field approximation is correct in the first order 
in membrane concentration n — S/2V for all frequencies, 
as justified in section on the origin of the y/co dispersion. 

The self-energy shifts the pole of the propagator 
thereby changing the diffusive dynamics due to disorder- 
averaged interaction with random membranes. Remark- 
ably, on this level there are no higher-order terms in pow- 
ers of q 2 . This allows us to represent interaction with 
membranes solely as renormalization of the effective dif- 
fusivity Do — > D (w), equation (4), which acquires fre- 
quency dependence dictated by that of the vertex (18). 
The role of 2?(w) is similar to that played by the disper- 
sive refraction index. Indeed, the relation (20) is analo- 
gous to the mean-field relation between the refraction in- 
dex and the forward scattering amplitude 40 . In our case, 
the presence of "scatterers" (membranes) introduces the 
memory kernel into the current-density response. 

In higher dimensions d > 1, averaging is performed 
both over positions of membrane in the direction n 
normal to its surface, and over its orientations. This 
amounts to a frequency shift iuj — » iw(q) = ico — -Doqf , 
£ w , q; „ = (Sy2T0V w(q) . qn , qn , with q|| = q (qn)n the 
conserved momentum parallel to membrane, and subse- 
quent orientational average S Uiq = (S w , q:n ) n - Substitut- 
ing the latter self-energy part in equation (19), obtain 
the result (4) for any d. Note that the y/u dispersion, 



a signature of the one-dimensional character of the scat- 
tering problem in the direction normal to the membrane, 
survives the orientational average in d > 1. 

Effective medium picture. The RG approach as- 
sumes that the boundary condition (1) is adjusted ac- 
cordingly, D — > T>{lo)\l at each RG scale L. Strictly 
speaki ng, this i s valid only for sufficiently small frequen- 
cies \JT>{uj) /oo 3> 1/n, corresponding to diffusing past 
many membranes, when the system can be treated as 
a uniform effective medium at each RG step. Luckily, 
for the opposite case of large frequencies the full result 
(5) matches the exact y/uFF >• ( limit (4), corresponding 
to the original boundary condition (1). This allows us to 
keep the frequency simply as a parameter in equation (8). 

From this discussion one expects that the agreement of 
the RG result (5) with the simulations is the worst in the 
intermediate regime, when y/vjuj ~ 1/n, equivalent to 
t <~ T r . Fig. 2a shows that the discrepancy between the 
RG result and the simulations is indeed highest around 
t <~ r r and remains under 25%. The agreement improves 
for t >• r r to within 10% even for the largest (. 

We note that in d = 1, the RG limit Ci(oo) overes- 
timates the e xact result for the one-dimensional lattice 
model 36-39 , \/2/tt, by the factor of two. We also observe 
this in the one-dimensional simulations (not shown). 
There are no exact results for this problem for d > 1. As 
it is often true with mean-field-like description, the ac- 
curacy of our solution is likely to improve for increasing 
d. This is consistent with the notably better agreement 
(Fig. 2c) of the theoretically obtained coefficient C*2(C) 
with that determined from the d = 2 simulations. 

Monte Carlo dynamics was realized on a square 
patch with periodic boundary conditions, embedded with 
randomly placed and oriented membranes, cf. Fig. 1, 
of nearly identical S/V. For the simulation of each £ 
in Fig. 2, the diffusion coefficient D{t) was calculated 
by averaging the displacement variance (a; 2 ), cf. equa- 
tion (3), over a total of 4 x 10 5 random walkers evenly 
split between 40 disorder realizations, with an average 
of 385 membranes per patch. The trajectory of each 
random walker was a sequence of moves in a randomly 
chosen direction over a distance dr = yj AD dt during a 
time step dt, with the total diffusion time t up to 100r r , 
corresponding to a maximum of 4 x 10 6 time steps per 
walker. Transmission across a membrane occurred with 
probability P cx Kdr/D -C 1, cf. ref. 26. The disorder 
strength ( was varied by changing k. The time step dt 
for each ( is chosen so that P < 0.007, and the ratio 
dr/a = Sdr/W < 0.1. We calibrated our results for 
a quasi-one-dimensional disorder (random parallel mem- 
branes), which reproduced the exact limit (2) with about 
1% accuracy. The random walk simulator was developed 
in C++. Simulations were performed on the NYU Gen- 
eral Cluster. With an average of 200 CPU cores used 
simultaneously, all simulations took about 100 hours. 
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